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Abstract 

We derive a series expansion for the multiparameter fractional 
Brownian motion. The derived expansion is proven to be rate op- 
timal. Keywords: fractional Brownian motion; series expansion; 
Bessel functions. 

1 Introduction 

The fractional Brownian motion with Hurst parameter H e (0, 1) is defined 
as the centred Gaussian process £(t) with the autocorrelation function 

R(s,t) = E£(8)t(t) = \{\s\ 2H + \t\™ -\s- t\' H ). 

This process was defined by Kolmogorov [3] and became a popular statistical 
model after the paper by Mandelbrot and van Ness [13] . 

There exist two multiparameter extensions of the fractional Brownian 
motion. Both extensions are centred Gaussian random fields on the space 
M. N . The multiparameter fractional Brownian sheet has the autocorrelation 
function 

1 N 

R (*>y) = ^nII^\ 2Hj + \y>\ 2Hj - \*j -vA 2Hj l E i e (M) 
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while the multiparameter fractional Brownian motion has the autocorrelation 
function 

i?(x )y ) = i(iixir+iiyir-iix- y m, ^ 

where || ■ || denotes the Euclidean norm in and where H e (0, 1). 

Dzhaparidze and van Zanten [3| and Igloi [7] derived two different explicit 
series expansions of the fractional Brownian motion. In [1] , Dzhaparidze and 
van Zanten extended their previous result to the case of the multiparameter 
fractional Brownian sheet. All the above mentioned expansions were proven 
to be rate optimal. We extend the results by Dzhaparidze and van Zanten 
to the case of the multiparameter fractional Brownian motion. 

To present and prove our result, we need to recall definitions of some 
special functions [5]. The Bessel function of the first kind of order v is 
defined by the following series 

J ^ = E 2 -m!i> + m + l)- (2) 

m=0 v ' 

Let < ju,2 < ■ ■ ■ < j u ,n < • • ■ be the positive zeros of J u (z). Let 

, m H = 2^- 2 >/ 2 r(Ay 2 ) ^ m>0, (3) 

and let 5^ denote the Kronecker's delta. 

The Gegenbauer polynomials of order m, C^, are given by the generating 
function 

v ' m=0 

Let m be a nonnegative integer, and let m , mi, . . . , m N _ 2 be integers 
satisfying the following condition 

m = m > rrii > • • > ^n-2 > 0. 
Let x = (xi, x 2 , ■ ■ ■ , xn) be a point in the space M. N . Let 



fk — \l ' x k+l + X k+2 + ' ' ' + x \i 
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where k — 0, 1, . . . , N — 2. Consider the following functions 

(, • \ ±mjv_2 N— 3 

XN -' + tXN rT 2 2 TT r T k ~ mk+1 



fc=0 

+(JV-fc-2)/2 / V 



and denote 

y(m fc , ±, x) = r m H(m k , ±, x). 

The functions Y(m k , ±, x) are called the (complex-valued) spherical harmon- 
ics. For a fixed m, there exist 

him N) = Vrn + N-2)(m + N-3)\ 

n{m ' > {N-2)\m\ [) 

spherical harmonics. They are orthogonal in the Hilbert space L 2 (S N ~ 1 ) of 
the square integrable functions on the unit sphere S 1 ^' 1 , and the square of 
the length of the vector Y(m k , ±, x) is 

, , yf 7i2 k - 2m *- N+2 T(m k ^ +m k + N-l-k) 

n ti ^ + ( N ~ 1 " fc )/ 2 )(^-i - ^) ! tr(m fe + (n - i - fc)/2)p 

Let Z = /(m fc , ±) be the number of the symbol (m , mi, ... , m N _ 2 , ±) in 
the lexicographic ordering. The real-valued spherical harmonics, S l m (x), can 
be defined as 

{Y(m k , +,x)/y/L(m k ), m N _ 2 = 0, 

\/2 Re Y(m k , +, x)/ y/L(m k ), m N _ 2 >0,l = l(m k ,+), 
-y/2lmY(m k , -,x)/ y/L(m k ), m N _ 2 > 0,1 = l(m k ,-). 

The hypergeometric function is defined by the series 



2 F 1 {a,b; c; z) = ^ 



-2 , 



(c) fc ife! 

where (u)^ = u(u + 1) . . . (u + k — 1), (m) = u. 
The incomplete beta function is defined as 

£ z (a,/3) = / t a -\l - tf- 1 di . 
Jo 
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Theorem 1. The multiparameter fractional Brownian motion £(x) has the 
form 



h(m,N) 



£( X ) = J2 Yl T ™n[g m (j\m-1\-H,n\\x\\) ~ S° m ]S l m ( m^m Hm ? 

m=0 n=l 1=1 ^ " " ' 

where ^ l mn are independent standard normal random variables, and 



(5) 



)//+! 



7r (iV-2)/2 r ( jf 7 + N/2)T(H + 1) sin(7r/J) 

r(^/2)j| m _ 1 |_ i?+ i(j| m _i|_ i?in H m + 1 ll _„ n 



(6) 



T/ie series ([5]) converges with probability 1 m t/ie space C(i3) o/ continuous 
functions in B = { x G : ||x|| < 1 }. 

Ayache and Linde p] derived another representation of the multiparam- 
eter fractional Brownian motion, using wavelets. Hence it is natural to com- 
pare their result with Theorem [U 

Consider the multiparameter fractional Brownian motion £(x), x G B 
as the centred Gaussian random variable £ on the Banach space C(B) of 
all continuous functions on B. The pth /-approximation number of £ [10] is 
defined by: 



l P (0 = mf 




3=P 



1/2 



C(B), 



i=i 



where £j are independent standard normal random variables and the infimum 
is taken over all possible series representations for £. 

Ayache and Linde [T] determined the convergence rate of l p (£) — > as 
p — > oo. To formulate their result, introduce the following notation. If a n , 
n > 1, and b n , n> 1 are sequences of positive real numbers, we write a n -< b n 
provided that a n < cb n for a certain c > and for any positive integer n. 
Then a n m b n means that a n ^ b n as well as b n ^ a n . In the same way, 
we write /(u) ^ fi'(u) provided that /(u) < cg(u) for a certain c > and 
uniformly for all u, and /(u) ~ #(u) if /(u) ^ g(u) as well as g'(u) ^ /(u). 

Ayache and Linde [U Theorem 1.1] proved that ~ p _ ^ 7V (logp) 1//2 . 
They also proved that their wavelet series representation possesses the opti- 
mal approximation rate. We will prove that our representation (jSJ) is optimal 
as well. 
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Theorem 2. The representation <^ possesses the optimal approximation 
rate for £ on B. 

Theorem [1] is proved in Section [2] while Theorem [2] is proved in Section [31 
I am grateful to Professors K. Dzhaparidze, M. Lifshits, and H. van Zan- 
ten for useful discussions. I thank the two anonymous referees for their very 
careful reading of the manuscript and for their helpful remarks. 



2 Proof of Theorem H 

It is well known p], that the multiparameter fractional Brownian motion can 
be represented as the stochastic integral 

e »(p,x) _ i 



r e HP> x J _ i 

£(x) = C H iV / v || n ||JV/2+ff dW (P 

Jr n IIPII ' 



where dW is the complex-valued white noise obtained by Fourier transfor- 
mation of the real-valued white noise. It follows that the covariance function 
P]) of the multiparameter fractional Brownian motion can be represented as 



i2(x,y) = c\ N I [e^ - l][ e -W> - l]||p||- M dp. (7) 

R N 



The following Lemma gives the explicit value of the constant c 2 HN . This 
result was announced by Malyarenko [T2] . 

Lemma 1. The constant c 2 HN has the following value: 



r 2 



2*B-ir(H + N/2)T(H + 1) sin(TriJ) 

^V+2)72 • 

Proof. For N = 1, our formula has the form 

2 2 2H - l T(H + l/2)r(# + l)sin(TriJ) 
c hi ~ ^ ' 

or 

2 T(2H + l)sin(TrF) 

C H1 — 7, • 

27T 

Here we used the doubling formula for gamma function. This result is known 
|17j . Therefore, in the rest of the proof we can and will suppose that N > 2. 



5 



Rewrite ([7j) as 

R(x,y) = c 2 HN [ [l-e^]||p|r^rfp 
Jn N 

+ c 2 HN [ [l-e-^Hpir^dp (8) 

-c 2 HN f [l - e^-^]\\ P \\- N -™ dp. 
Jn N 

Consider the first term in the right hand side of this formula. Using formula 
3.3.2.3 from HU, we obtain 



c 2 HN [l-e^llpll-^dp 
Jr n 



HN 



POO P7T 

/ X N ~ 1 d\ / [l-e lX ^ cosu ]X- N - 2H sm N - 2 udu. 
Jo Jo 



r((iv-i)/2) J 

It is clear that the integral of the imaginary part is equal to 0. The 
integral of the real part may be rewritten as 



r 2 



I [l-e^Hpll-^-^dp 
Jn N 

o (jV-l)/2 2 roo ptt 

= r((JV - m) I l^-^Mcosu)]sin N - 2 ududX. (9) 

To calculate the inner integral, we use formulas 2.5.3.1 and 2.5.55.7 from 



' . N _ 2 V^T((N - l)/2) 
8111 UdU = fp) ' 

cos(A||x|| cos M ) sin^" 2 ^ = ^ N ^I 2 T{{N - l)/2) ^Taw • 
o (^ll x ll) ( )7 

It follows that 

^[l-cos(A||x||cos M )]sin jV - 2 M ^= ^^"^^ [l-go(A||x||)]. (10) 
Substituting ffTOl) in (jUJ), we obtain 

c^y R Jl-e^]||p||-^d P = ^^ jf A- 2 ^[l-^o(A||x||)]dA. 

(11) 
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To calculate this integral, we use formula 2.2.3.1 from 

1 (1 _ t , 2) ,«- 3 ,/ 2(itl= Ar(( J v-i)/2) 



It follows that 



2T(N/2) 



On the other hand, according to formula 2.5.6.1 from [T3] we have 



\l - v*r-^ cos(A||x||„) dv = ^2^F((N - l)/2) ^ffji } 



It follows that 



^o(A||x||) = r-^!! l2 \ /0 , / (1 - v^ N ^ cos(\\\x\\v) dv. (13) 



y/*T{(N-l)/2) Jo 

Subtracting f|T3|) from ([12"]) . we obtain 

l-gb(A||x||)= ^Z { H ,2 \^ f '(I - ^ 2 ) M/2 sin 2 



o 



0Fr((iv - 1)/2) 

Substitute this formula in (ITT]). We have 



Jr n 



-(N-l)/2 2 /-oo /•! 



r((iv-i)/2)7 y 

Consider two integrals: 



A -2H-l (l _ u2) (iV-3)/2 sin 2 ( ) dvdX (M) 



1 A- 2 ^- 1 sin 2 ( *^a! dA and /°° A" 2 "" 1 sin 2 ( ^A ) <1\. 



'0 \ * / ./l 

In the first integral, we bound the second multiplier by A 2 /4. In the second 
integral, we bound it by 1. It follows that the integral 

oo / II I 

X -2H-1 gin 2 / «||X|| A I ^ 



converges uniformly, and we are allowed to change the order of integration in 
the right hand side of (1141) . After that the inner integral is calculated using 
formula 2.5.3.13 from \l l 



o 



r 3,M snrf^l ^A = - ™™ , rr J |x" 3 " 
2 J 4T(2H + l)sin(Tr^) 



Now formula ( TTTT) may be rewritten as 



c 2 HN [l-e^)]||p||«rfp 
Jn N 



I v ^^_ v 2 ){ N-3)/2 dv „„||2H 



X 



T((N - \)/2)V{2H + 1) sin(7r/J) J 
The integral in the right hand side can be calculated by formula 2.2.4.8 from 



2H 



V 



r((Jv-i)/2)r(ir + i/2) 







2T(H + N/2) 
and (TTTJ) is rewritten once more as 



c 2 HN [l-e^llpll-^dp 
Jr, n 



\2H 



(iV+2)/2 2 

= 71 C HN ||xll 2H (^^) 

2 2H T {H + 1)T {H + N/ 2) sin{7iH) 11 11 ' y } 

On the other hand, the left hand side of (fT5l) is clearly equal to — 1| j 
The statement of the Lemma follows. □ 

Taking into account (fTTj) . one can rewrite (jSJ) as 

r, JV/2 2 poo 

i?(x,y) = — — ^ / A- 1 -^[l-^o(A||x||)-^ (A|| y ||)+^o(A||x-y||)]rfA. 



T(iV/2) J 

(16) 

Let x ^ 0, y ^ and let <p denote the angle between the vectors x and y. 
Two addition theorems for Bessel functions (formulas 7.15(30) and 7.15(31) 
from [61 vol. 2]) may be written in our notation as 



oo „(N-2)/2, \ 

g (X\\x - y||) = Km, iV)g m (A||x|l) gm (A||y||) m ^ J> , (17) 

dm ( 1 J 



m=0 



where g m is as in fl3j). Substituting (fT7|) in (ITH|) . we obtain 

m=0 

where 



9 n/2 o° r (7V ~ 2)/2 iWi/^ z 100 

^ (X ' y) = E fe(m ' ^ 7 «WlW^||(A)dA, (18) 



C(A) = CHivA- H - 1/2 [0 m (rA) - <]. 

Recall that the Hankel transform of order v > — 1 of a function k G 
L 2 (0, oo) (see [HI Section 8.4]) is defined as 

POO 

k(u) = / k(X) J u {uX)\fuX dX. 
Jo 

Define the kernel k™(u) as the Hankel transform of order \m — 1| — H of the 
kernel k™(A): 



k?(u)= K™(X)J lm - lhH (uX)VuXdX. (19) 

Therefore, the Parceval identity for the Hankel transform [181 Section 8.5, 
Theorem 129] implies that (fT8|) may be rewritten as 

N /2 °° r (iV-2)/2 r x 

%-^ W E^p^H^)| ^(t,)^,,^)^ (20) 
To calculate £™(u) in the case of m > 1, we use formula 2.12.31.1 from 

"V W ~ 2 H+N / 2 - 1 T(H + N/2)r m+N - 2 X( -°' r ^ u) > ^ > 

where X(o,r)( u ) denote the indicator function of the interval (0,r). 
For m — 0, the integral in (TT9"j) can be rewritten as 



k° r (u) = 2^l 2 T{N/2)c HN r^ N l 2 ^ / A 1 "^ 2 J ( ^_ 2)/2 (rA) Jx-h{uX) dX 

Jo 

POO 

- y/uc HN / X~ H Ji_ H (uX)dX. 
Jo 

To calculate the second integral, we use formula 2.12.2.2 from [T5] . 
[^ H J^n{nX)dX= 



For the first integral, we use formula 2.12.31.1 from [15J once more. In the 
case of u > r we obtain 

rn - mA N - 2 )/ 2 



2H+{N-2)/2 T ( N l2) u l-H- 

In the case of u < r, we have 

T(l - H)r 2H+N/2 ~ 3 



\ l - H - N l 2 J^ N _ 2 y 2 (r\)J 1 _ H (u\) dX 



2 h+(n-2)/2 V ( H + (jv - 2)/2)r(2 - H) 
x 2 F 1 (l-H,2-H- N/2; 2-H; u 2 /r 2 ). 

Using formula 7.3.1.28 from |16] , we can rewrite the last expression as 

r (N~2)/2 

X l-H-N 2 J (JV _ 2)/2 ( rA ) J^uX) dX = 2 H + (N-2)/2 T ( H+ ( N _2)/2) U l-V 

x B u 2 /r2 (l- H,H + (N -2)/2). 
Combining everything together, we obtain 

x B u2/r2 (l -H,H + {N- 2)/2)]x ( o,r)(w). 

It follows from the last display and from (j2ip that the support of the kernels 
^fxliC") anc ^ ^llyll^) ^ es m s i nce ll x ll) l|y|| < 1 for x, y G B. We rewrite 
(1201) as 

n N/2 °° nW-Wfcrnm} f 1 
\ ' ' m=0 ( - ym 1/0 

for x, y G £>. 

For any v > — 1, the Fourier-Bessel functions 



= 1 j^—rMj^nU), n > 1 

form a complete, orthonormal system in L 2 [0,1] [201 Section 18.24]. Put 
v = \m — 11 — H. We obtain 



)du = J2K(MK l (\\y\\), 



n=l 
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where 

-i 



K(M)= / K%\\{u)W\m-l\-H,n{u) du. 
JO 

This integral can be rewritten as 

POO 

C(NI)= / *T*\\(u)V\rn-l\-HMdu. 
JO 

To calculate it, we use the inversion formula for Hankel transform: 

k™{u) J v {Xu)Vhldu = K?(\) 



and obtain 



v/2 



C(l|x||) = — . ^OlM-^)- 

J|m-l|-i?+lU|m-l|--ff,nJ A/J|m-l|- " 



Substituting the calculated values in (1221) . we obtain 

y) - FW2) E E ,!(m ' iV) 

K lixll U|m-l|— ^nj^lkrll U|m-l|— i?,n) 



(23) 



^|ra-l|-jf+l(i]m-l]-#,rc)j|m-l|-If,ri 



According to the addition theorem for spherical harmonics [6, Vol. 2, 
Chapter XI, section 4, Theorem 4], 

{N-2)/2, s h{m,N) 



£*(h)*(h)- (24) 



C^- 2)/2 (l) " T(N/2)h(m,N) £j 
Substituting this equality in ( 1231) . we obtain 

_ iu oo oo h(m,N) m /■ \ m /■ \ 

(x ' yj "(r(iv/2))^ o ^ 



iy 
n 



It follows immediately that the random field £(x) itself has the form ([5]), (151), 
and the series (JSJ) converges in mean square for any fixed x 6 B. 

Since the functions in ([5]) are continuous and the random variables are 
symmetric and independent, the Ito-Nisio theorem [19J implies that for prov- 
ing that the series (jSJ) converges uniformly on B with probability 1, it is 
sufficient to show that the corresponding sequence of partial sums 

M M h(m,N) , v 

£m( X ) = 5^S Yl T mn[g m (j\ m—l\ — H,n\\~X-\\^) ^m\^m ( i i ) ^mn (^^) 
m =0 n=l 1=1 ^" 

is weakly relatively compact in the space C{B). To prove this, one can use 
the same method as in the proof of Theorem 4.1 in [5]. Proof of Theorem [1] 
is finished. 

Consider some particular cases of Theorem [TJ In the case of N = 1, 
there exists only two spherical harmonics on the O-dimensional sphere S° = 
{— 1, 1}. They are 

So( x ) = ~^=j Sii x ) = ~/^> 
where x G S°. Moreover, we have 

9o(jl-H,n\x\) = COs(ji_H,nM), 9l{j-H,n\M) = ^{j-H,n\x\). 

Therefore formula (jSJ) becomes 

On ( f cos(ji- g , n x) - 1 x ^ S in(j- g , n x) x 

\ n=1 J2-H[3l-H,n)3i-H,n n= \ J l-H{3-H,n)3-H,n 

(26) 

This is the result of [3]. 

Consider the case of the multiparameter fractional Brownian motion on 
the plane (N = 2). Let (r,(p) be the polar coordinates. The spherical 
harmonics are: 

SJM = 4=. S» = ^, S»- Sin(m<rf 



2' lib \ I / / ' //£• \ < / / 

7T V7T V 71 " 
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It follows that 



2 H + l r(H + l)^sm(nH) ( 1 ~ J (ji-H, n r) - 1 x 



oo oo 



\ - \ - J m (jm-i-H, n r) cos(my) x 

m=1 n=1 ■Jm-H\Jm-l-H,n)J m -i-H,n 

J m (j m -i-H, n r) sin(myg) 

m= l n=1 Jm-H\Jm-l-H,n)J m -i-H 



2 

van 



3 Proof of Theorem [2 



It is enough to prove that the rate of convergence in §5§ is not more than 
the optimal rate p~ H l N (\ogp) l l 2 , where p denote the number of terms in a 
suitable truncation of the series. Since for N > 1 we have a triple sum in our 
expansion (jSJ), it is not clear a priori how we should truncate the series. We 
need a lemma. 

Lemma 2. The Bess el functions, the functions g m defined by (j3J), the num- 
bers h(m,N), and the spherical harmonics ^(x/HxH) have the following 
properties. 

1- 3u,n ~n + v/2- 1/4. 

3. \g {u)\ < 1. 

4- \9m(u)\ d [m(m+JV _ 1 2)] (iv- 1 )/4 , m > 1, JV > 2. 

5- |<4(u)| d [m(Tn+J v- 1 2)] (iv-3)/4 , m > 1, N > 2. 
5. h(m,N) ^ m^" 2 , n > 2. 

7. |^(x/||x||)| ^m^" 2 )/ 2 , N > 2. 

Proof. Property [H is proved in [201 Section 15.53]. It is shown in [20J, Sec- 
tion 7.21] that 

J 2 M + Ju+i («)«-, (27) 
w 
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since Property [21 Property [3] follows from the fact that go(u) is the element 
of the unitary matrix |21j . 

Let \i\ < //2 < ■ ■ ■ be the sequence of positive stationary values of 
the Bessel function J m +(N-2)/2(u). It is shown in [20j Section 15.31] that 

\Jm+(N -2) /2(^l)\ > \J m +(N — 2)/2 (^2) | > 

Let Mi be the first positive maximum of the function g m {u). Let x\ denote 
the maximal value of the function \g m {u) \ in the interval (0, j m +{N-2)/2,i), let 
X2 denote the maximal value of |<7 TO («)| in the interval (j m +(N-2)/2,i, jm+(N-2)/2,2), 
and so on. Then we have 

\Jm+(N -2) /2 (/^l) I 



x 2 d> 
X3 d 



(iV-2)/2 
1 



Jm+(N -2) /2(f J '2) I 
.{N-2)/2 
^m+(AT-2)/2,l 

^m+(V-2)/2(W3)| 
.{N-2)/2 
3 m +(N-2)/2,2 



and so on. The right hand sides form a decreasing sequence. Then 

\Jm+{N -2) /2 vA*l) I 



(JV-2)/2 
1 



It follows from (E7D that 



and we obtain 



|J„(«)|=< ' 
<u 



\g m (u) \ d 



Wl ^- 2)/2 



To estimate consider the differential equation 
f (if 

m 2 ^4 + {N- l)u-r- + b 2 - m(m + iV — 2)1/ = 0. (28) 
aw 2 aw 

It follows from formulas (3) and (4) in [201 Section 4.31] that the function 
g m satisfies this equation. It follows from the series representation (J2J) of the 
Bessel function and the definition ([3]) of the function g m (u) that for m > 1 
we have g m (0) = and g m increases in some right neighbourhood of zero. 
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Therefore we have g m {u\) > 0, g' m {ui) = and g'^iui) < 0. It follows from 
equation (|2"B"|) that U\ > y/m(rn + N — 2). 

For Hi, we have the estimate \i\ > m + (N — 2)/2 ([20, Section 15.3]). 
Using the inequality m+(N— 2)/2 > yTnjm + iV — 2), we obtain PropertyHl 

In any extremum of g' m {u) we have g'^iu) = 0. It follows from (1281) that 



[m 2 — m(m + iV — 2)]p m («) 



(JV - l)u 



Property [5] follows from this formula and Property HI 
Property [6] follows from (jlj) and Stirling's formula. 
It follows from §M§ that 



S l 



2 < T(N/2)h(m,N) C { I- 2)/2 (cos<p) 



2tt n / 2 



a 



(N-2)/2, 



By P Vol. 2, formula 10.18(7)], 



max |C(f- 2 )/ 2 (cos^)|=C(f- 2 )/ 2 (l). 



It follows that |S^(x/||x||)| ^ ^//i(m, iV). Now Property [7] follows from 



Property [6l 
Denote 



□ 



«L( X ) = 7"mn|>m0'|m-l|-ff,n||x||) - <^JS^ 



Using Lemma [2, we can write the following estimate: 

-i 2 



h(m,N) 

^ ] U mn( X )C 
1=1 



I 

mn 



1 



(m + l)(m/2 + n) 



2H+1 



For any positive integer q, consider the truncation 

h(m,N) 

(rn+l)(m/2+n) 2H + 1 <q 1=1 



(29) 



(30) 
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of the expansion (151) . The number of terms in this sum is asymptotically 
equal to the integral 



u du dv 

(u+l)(u/2+v)<q,u>0,v>l 

and therefore is bounded by a constant times q N /( 2H+2 \ Therefore we need 
to prove that 



/ 



h(m,N) 

E E 

(m+l)(m/2+n) 2H + 1 >q 1=1 



\ 



C(B)J 



1/2 



(N/(2H+2)){~H/N) 



(\ogq 



AT/(2if+2)\l/2 ^ y-H/(2H+2) 



'(logg) 



1/2 



By equivalence of moments ([9j Proposition 2.1]), the last formula is equiva- 
lent to the following asymptotic relation. 



Esup 



(m- 



E 

-l)(m/2+n) 2H + 1 



>q 



h(m,N) 

E 

1=1 



-< 



-H/(2H+2) 



(log q) 



1/2 



To prove this relation, consider the partial sum ?7a;(x) defined by 

h(m,N) 

%( x ) = ^ ] ^ ] U mni y ^)^,mn- 

2 k - 1 <(m+l)(m/2+n) 2H + 1 <2 k 1=1 



For a given > 0, let x 1; . . . , xp fc e 5 be a maximal e^-net in S , i.e., 
the angle 99 (x^, x&) between any two different vectors Xj and x& is greater than 
Ek and the addition of any new point breaks this property. Then Pk ~ 1 . 
A proof of this fact for the case of N = 3 by Baldi et al [21 Lemma 5] , is easy 
generalised to higher dimensions. The Voronoi cell S(~x.j) is defined as 

Sfe) = { x G S 1 *' 1 : <p(x,Xj) < v?(x,x fc ), A; }. 

Divide the ball B onto [e^ 1 ] concentric spherical layers of thickness ^ 
Voronoi cells determine the division of each layer onto P^ segments. The 
angle y?(x, y) between any two vectors x and y in the same segment is ^ 
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(we never choose the point 0). Call all segments in all layers Bj, 1 < j < 
Mf, = -Pfc^fc 1 ] di £ k N ■ Choose a point Xj inside each segment. Then we have 

Esup |?7fc(x)| < E sup \r) k (xj)\+E sup sup |%(x) - r] k (y)\. 

xeB l<j<M k l<j<M h x.yeBj 

By a maximal inequality for Gaussian sequences (22J Lemma 2.2.2], the 
first term is bounded by a positive constant times 



V'l + logMfc sup JE(ji k (xj))' 
i<j<M k v 



Using the estimate (1301) . we obtain 



E(^(x)) 2 z< If — — 

J J2 k - 1 <(u+l)(u/2+v) 2H + 1 <2 k ,u>0,v>l \U-\- I 



du dv 



)(u/2 + v 



\2H+1 ■ 



The integral in the right hand side is easily seen to be ^ 2 kH /( H + 1 ) m it 
follows that 



E sup |%(x,.)| <2- fe ^ +2 )v/l + logM fc . (31) 

l<j<M k 

To estimate the second term we write 
E sup sup \rj k (x) - rj k (y)\ 

l<j<M k x,yGBj 

h(m,N) 

<E SUp SUp Yl l M mn( X ) - «L(y)l • I^Ll- 

l<j<M k x.yeB, 2 fc-i <(m+1)(m/2+n) 2 H+ i< 2 fc i=1 

The difference \u l mn {x) — u l mn (y)\ can be estimated as 

l«L(x) " «L(y)l < \9m(J\m-l\-H, n \\ X \\) - ■ |^(x/||x||) - S l Jy / \\y \\) \ 

+ i^(y/l|y||)i-IUii 

m-l|-Jf",n||x||) — g m (j\ m—l\—H,n || y || , 

Let A denote the angular part of the Laplace operator in the space W N . 
Using the Mean Value Theorem for g m , Lemma [2] and formulas 

AX(x/||x||) = -m(m + N- 2)^(x/||x||), 
|^(x/||x||) - S* m (y/||y||)| d sup |/=A^ l (z/||z||)| ¥ >(x,y), 
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where the sup is taken over all points z belonging to the segment of the 
geodesic circle connecting the points x/||x|| and y/||y||, we obtain 

l«L( x ) - «L(y)l =< e k T mnV ^{m/2 + n). 



Hence, we have 



E ^sup^ sup^ |77 fc (x)— »y fc (y)| ^ e fc h(m,N)T mn ^/m(m/2+n). 

2 k - 1 <{m+l){m/2+n) 2H + 1 <2 k 



l<i<M fc x,yeBj 



Using Lemma [2] once more, the last formula may be rewritten as 



E_sup^ supjr/ fc (x)-^(y)| < e k ^ m iV - 3/2 (m/2+n) 1/2 - ff . 

2 k ~ 1 <(m+l)(rn/2+n) 2H + 1 <2 k 



l<i<M fe X ,ye&/ 



The integral, which corresponds to the sum in the right hand side, is easily 
seen to be asymptotically equal to 2 kN ^ 2H+2 \ Therefore, 



E sup sup |%(x) - Vk (y)\ * e k ■ 2 kN ^ 2H+2 \ 

l<i<M fc x,yeBj 



Combining this with the estimate fl3Tl) for the first term, we obtain the asymp- 
totic relation 

Esup | % (x)| ± 2- kH ^ 2H+2 ^l + \ogM k + e k ■ 2 kN ^ 2H+2 \ 

NOW Set S k = 2-HH+N)/(2H + 2) and recaU th&t Mfc ^ £ -N = 2 kN(H+N)/(2H+2) _ 

Then we see that the first term is bounded by a constant times 2~ kH ^ 2H+2 ^\/k 
and the second one is of lower order. We proved that 

Esup |%(x)| z< 2~ kH ^ 2H+2 ^Vk. (32) 



To complete the proof, it suffices to show that 



E sup 

xGB 



h(m,N) 

E E 

(m+l)(m/2+n) 2H + 1 ><j 



'It.' 



(x)d 



-H/(2H+2) 



(logg) 



1/2 



(33) 
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Let r be the positive integer such that 2 r 1 < q < 2 r . Then by the triangle 
inequality 



iinn 



< 



h(m.N) 

E E we. 

(m+l)(m/2+n) 2H + 1 >g i=l 

h(m,N) 

U rnn( X -)^mn 

q<{m+l)(m/2+n) 2H + 1 <2 r 1=1 

h(m,N) 

E E 

g<(m+l)(m/2+n) 2H + 1 <2 r J=l 



h(m,N) 

E E we. 

(m+l)(m/2+n) 2H + 1 >2 r 1=1 



mn 



<X>*wi + 



k>r 



Using (|32[) and the fact that 



k>r 



we obtain 



^Esup|7? fc (x)| ^ ^v^2- 



■kN/(2H+2) 



k>r 



k>r 



^ g -if/(2ff+2) (log?) l/2_ 

The arguments in the proof of fl32l) show that since 2 r ~ 1 < g < 2 r , we also 
have 



E sup 



h(m,N) 

/ . / , " 1 mn\ ftmn 

g<(m+l)(m/2+n) 2£r + 1 <2'' Z=l 



^ g -^/(2H+2) (logg) l/2 



Relation (1331) is proved. 
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